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Abstract 

We numerically investigate the jet propagation through a rotating collapsing Wolf-Rayet star with 
detailed central engine physics constructed based on the neutrino-driven coUapsar model. The col- 
lapsing star determines the evolution of mass accretion rate, black hole mass and spin, all of which 
are important ingredients for determining the jet luminosity. We reveal that neutrino-driven jets in 
rapidly spinning Wolf-Rayet stars are capable of breaking out from the stellar envelope, while those 
propagating in slower rotating progenitors fail to jet breakout due to insufheient kinetic power. For 
progenitor models with successful jet breakouts, the kinetic energy accumulated in the cocoon could 
be as large as ^ lO^-'^erg and might significantly contribute to the luminosity of the afterglow emis- 
sion or to the kinetic energy of the accompanying supernova if nickel production takes place. We 
further analyze the post breakout phase using a simple analytical prescription and conclude that the 
relativistic jet component could produce events with an isotropic-luminosity Lp(^iso) ^ lO^^erg/s and 
isotropic-energy Eji^i^g-^ ^ lO^^erg. Our findings support the idea of rapidly rotating Wolf-Rayet stars 
as plausible progenitors of GRBs, while slowly rotational ones could be responsible for low luminosity 
or failed GRBs. 

Subject headings: supernovae: general — neutrinos — hydrodynamics, black hole physics 



1. INTRODUCTION 

Long-duration Gamma-Ray Bursts (GRBs) are 
thought to originate fro m the death of massive stars 
(jWooslev fc Blooml |2006[ ). It is widely recognized that 
the study of GRBs provides the important knowledge 
on the final evolutionary stage in the life of massive 
stars. Although the nature of GRBs remains elusive, one 
viable scenario to produce a GRB is the neutrino-driven 
collap sar model (jWooslevI 119931 : iMacFadven fc Woosleyl 
Il999( ) . The gravitational collapse of the rapidly rotating 
core is believed to create a fast rotational Kerr black 
hole. The copious amounts of neutrinos and their anti 
particles, which are emitted from the hot accretion 
disk, annihilate and create a n electron positron fireball 
around the rotational axis (lEichler et all [l989l) . The 
baryon-starved fireball is believed to give rise to a 
relativistic coUimated outflow, and eventually produce 
a GRB after the beam has broken free from the stellar 
progenitor. 

Over the years, substantial work has been made 
towards understanding if neutrino-driven collapsar 
jets c an produce the required rela. ti yistic outflow (see 
e.g. lMa.cFadven fc WooslevI (119991): iFryer & Meszaros 



central engine is roughly six orders of magnitude longer 
than the dynamical timescale of the nascent black hole. 
Computational resources are not yet available to perform 
such numerical studies. Because of these difiiculties, a 
number of studies have thus far emplo ye d steady state 
approximations JPopham et al.' 1999: 'Naravan & Yjl 
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iTavlor et all (|20Hl Vr However, there are still 
many open questions. In terms of the central engine, 
one of the largest uncertainties is the efficiency of energy 
deposition by neutrinos. Although numerical studies 
are very powerful methods to investigate the energy 
deposition rate by neutrinos and subsequent evolution of 
the jet, we need to solve General Relativistic Neutrino 
Radiation Hydrodynamics with microphysics for several 
tens of seconds after the black hole formation. This 
is certainly challenging as the typical life time of the 



Janiuk et ail 120071: IShibata et all l2007l : iHarikae et alj 
20101) with simplified microp hysical treatments. It 
is int eresting to note that IZalamea fc Beloborodovl 
()2011l) recently have conducted General Relativistic Ray 
Tracing Neutrino Radiation Transfer and they found 
that the energy deposition by neutrinos could be well 
described as a simple analytic formula. Thanks to these 
studies, the jet luminosity can be estimated qualitatively 
without the need for expensive numerical simulations. 

It should be noted, however, even if the neutrino- 
driven jet is successfully launched, this docs not guaran- 
tee the production of a GRB. As a minimum requirement, 
the jet needs to successfully penetrate the stellar enve- 
lope otherwise it would become non-relativistic and thus 
incapable of producing GRB. Ever since the neutrino- 
driven collapsar model was proposed, a number of nu- 
merical wor ks of jet propagation have been done (see e.g . 
(lAlov et al.ll20q0l:IZhang et al.l2003HMorsony et al.ll2007l: 
iMizuta fc Alovl 120091: iNagakura et al.ll2011l) ). IT these 
studies, however, the jet luminosity was assumed for sim- 
plicity to either be constant or to directly fol l ow th e 
mass accretion r ate (s ee e.g. IMacFadven et al.l ()200lD : 
INagakura et al.l (|2012[ )). In order to judge whether 
neutrino-driven jets are capable of successfully breaking 
out of their progenitors, and to explore the effects of ro- 
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Fig. 1.— SAM distribution of models (Red: Mref, Green: M150, 
Blue: M70, Pink: M50). The cross marks denote the intersec- 
tion between the SAM of star and SAM at the last stable orbit 
for a black hole with the mass and angular momentum inside the 
indicated coordinate. 

tation, one needs to take into account how the jet power 
scales with the neutrino energy deposition generated by 
the accompanying accretion. In this paper, we present, 
for the first time, the propagation of neutrino-driven jets 
employing ac curate neutrino energy deposit ion rate as 
calculated by iZalamea &: Beloborodovl ()2011[ ) . The evo- 
lution of mass accretion rate and black hole mass and 
spin, all of which are necessary to evaluate the energy 
deposition by neutrinos, are evaluated here using an in- 
ner boundary condition in the simulation, although the 
current calculation is still not self-consistent as it fails to 
resolve the accretion disk and assumes no feedback (e.g. 
a disk wind). The purpose of this study is (1) to clarify 
whether neutrino-driven jets can successfully break out 
from the stellar surface and (2) to determine the progen- 
itor's rotation rate necessary for successful jet breakout. 
As we will show in this paper, the final outcome of the 
explosion depends sensitive on the rate of rotation, and 
differences in rotation rate could be responsible for the 
observational differences seen between GRBs, low lumi- 
nosity GRBs (LLGRBs) and failed GRBs. 

2. METHODS AND MODELS 

We perform two dimensional, relativistic hydrodynam- 
ical axisymmetric (and also equatorial symmetric) sim- 
ulations of the accretion and subsequent jet propaga- 
tion. The numerical code employed in this paper is 
essentially the same as tho se used in previous papers 
(jNagakura et al.ll201lL l2012l ). The ini tial stellar density 
distrib ution is fixed as 16TI model in IWooslev fc Hegerl 
(|2006[ ). As is the case with previous studies, we cut the 
inner portions of the star from a certain radius. The self- 
gravity of matter in the active numerical regions is cal- 
culated by solving Poisson equations and the monopole 
gravity is added as the point mass at the inner excised 
region. The mass accretion rate (Af ) is estimated by 
th e mass flows thr o ugh t he inner boundary (see Eq. (1) 
in iNagakura et ahl ()2012[ )V The mass and angular mo- 



mentum in the excised region are assumed to be the same 
mass and angular momentum of the black hole. The time 
evolution of mass and spin of black hole is calculated by 
integrating mass and angular momentum flux crossing 
the inner boundary. It should be noted that when the 
specific angular momentum (SAM) at the location of in- 
ner boundary in equatorial region becomes larger than 
SAM at the inner most stable circular orbit (ISCO) (see 
cross marks in Figure [1]), we alter our prescription when 
calculating the angular momentum to: 



J a = M X J/sco, 



(1) 



where /q and Jisco denote the angular momentum flux 
and SAM at the ISCO, respectively. This treatment 
comes from the fact that the angular momentum of mat- 
ter in the disk is transported outwards due to the turbu- 
lent viscosity or non-axisymmetric waves, and finally the 
matter falls int o a black hole with ~ Jisco- 

According to IZalamea fc Beloborodovl (|201lh , the jet 
luminosity by the neutrino process is determined by; 



L, = 1.1 X 1052:,-4.8( i||^)-i 

^ (M < M,<,„) 

Xlrh-i {Ahgn < M < Mtrap) 

>Lp {M> Mtrap) 



(2) 



where rh = M/{Mq/s), Xms = rms/C^GMbh/c^) {vms 
denotes the marginally stable orbit). G and c de- 
note the gravitational constant and the speed of light, 
respectively. The characteristic mass accretion rate 
Mign and Mtrap are given as a function of the vis- 
cous parameter (a), an d the Kerr parameter (a) (see 
IZalamea fc Beloborodovl (|2011h ). Throughout this pa- 
per, we set a = 0.1 and Kerr parameter dependence on 
these accretion rate is linearly interpolated by Xms- 

The spherical symmetric density distribution is 
mapped by the spherical coordinate. The computational 
domain covers from lO'^ to 4 x 10^"cm. Note that the 
location of the inner boundary in the present study is 
ten times sm aller than in previou s jet propagation- stud - 
ies (see e.g. iLazzati et al.l (|2009D : iMizuta et alj ()201lD : 
INagakura et all (1201111 )7 The evolution of mass accretion 
rate , which is sensitive to t he location of inner boundary 
(see INagakura et all ()2012[ )). can thus be better captured 
by our simulations. However, as a result, these simula- 
tions become rather computationally expensive and we 
are only able to conduct them until the jet bow shock 
reaches the stellar surface or the black hole mass reaches 
1OM0. The evolution of the the post-breakout phase is 
then analyzed by using a simple analytic formalism (see 
Eqs. ©-(O). The results of an extended numerical sim- 
ulation will be nonetheless compare with the analytical 
approach for the reference model in order to confirm that 
the analytical approach qualitatively captures the evolu- 
tion of the jet dynamics (see Section IX^ and Figure Sj). 

We employ the gamma-law equation of state with 7 = 
4/3. The jet injection parameters such as the Lorentz 
factor and the specific internal energy are the same as 
those used in the standard model of INagakura et al.l 
()2012[ ). where the initial Lorentz factor and specific in- 
ternal energy are fixed to F = 400 and e = 0.01, respec- 
tively. It should be noted for a fixed F a nd e, the overall 
jet dynamics depend solely on dop (see INagakura et al.l 
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TABLE 1 

Summary of our results 



Model 


Breakout 


u 


^br 


Lip 








-^'j>f49.5 




^i>i-49..'i 








(s) 


(10^° erg/s) 


(lO^ierg) 


(10^1 erg) 


(10^1 erg) 


(lO^^erg) 




(s) 


Mref 


yes 


10.9 


27.8 


1.9 


1.4 


7.4 


4.4 


5.6 


27.8 


46.2 


M150 


yes 


2.2 


17.6 


3.2 


1.6 


11.7 


8.8 


9.8 


40.0 


55.0 


M70 


yes 


15.5 


45.5 


1.0 


0.9 


3.6 




1.7 




21.7 


M50 


no 


21.6 




0.6 















(|2012t l). In this study, we assume Oop — 9°, which agrees 
well with the opening angles deduced for Long GRBs 
(|Goldstein et al.ll2011[ ). The dependence of our results 
on the Oop will be discussed in section |3l 

The 1000 non-uniform radial grids cover all the compu- 
tational region while the meridian section is covered by 
60 uniform grids. The 3 level Adaptive Mesh Refinement 
techn ique, similar to that used in Nagakura_et al. (2011, 
l2012t l is also employed in order to decrease computational 
cost. We set up the stellar rotation by a sim ilar manner 
as those used in iLopez-Camara et al.l ()2010f ) . The SAM 
distribution is separated into radial and polar compo- 
nents as J(r, 6) — j{r)Q{6), where r and 6 are the spher- 
ical radius and polar angle, respectively. In the reference 
model (Mref), j(r) is given by 16TI model. For M150, 
M70 and M50 models, j(r) is multiplied by 1.5, 0.7 and 
0.5, respectively (see Figure [T] for the SAM distribution 
of our models) . The polar angle components are assumed 
to be rigid body rotation on shells, i.e., 0(6*) = sin^ 9. 

It is important to highlight that the simulations in 
this study do not cover the black hole accretion disk 
system. Even if the simulations cover the full com- 
putational domain, our numerical calculations can not 
treat the disk evolution appropriately, since the gen- 
eral relativistic effect and microphysics, which are im- 
portant in determining the disk evolution, are not in- 
co rporated. However, the analyt ical formula proposed 
by IZalamea fc Beloborodovl (|2011[ ) allows us to estimate 
neutrino luminosity without resolving the black hole ac- 
cretion disk system. Owing to these prescription, we can 
study the jet penetration phase as determined by the 
neutrino-driven energy injection. The study of the cou- 
pling between black hole and disk accretion system is the 
beyond the scope of this paper. 

We also note that the injection of the jet is de- 
layed in the simulation as a result of the inner core 
not poss essing eno ugh angular momentum to create 
a di s k (IWoosle-^fc Heger 20061: iLee fc Ramirez-Ruizj 
[200l IZalamea fc Beloborodovll2011f ). Based on the 
standard neutrino-driven coUapsar model (and as - 
sumptions those used in Chen fc Beloborodovl ()2007t ): 
IZalamea fc Beloborodo 

3 (I2011D '|. the central engine 
starts to operate after the accretion disk is formed around 
a black hole. Therefore, in our simulations, the jet is in- 
jected only when the SAM of matter at the inner bound- 
ary in equatorial plane exceeds Jisco- 

3. RESULTS 

3.1. Basic features and jet penetrability 

The overall evolution of the collapse of the progeni- 
tor seen in our simulations is not surprising ly similar to 
that found in iNagakura et al.l ()2011l 120121 ). The infall 
of stellar envelope generates a rarefaction wave, which 
propagates outwards. The envelope contraction is almost 



M150 



Mref 



M70 





Fig. 2. — The density contour (log scale) in the meridian section 
for each model at t = 17.6s (left), t = 27.8s (middle) and t = 46.5s 
(right). Each time corresponds to the time of jet breakout for 
M150, Mref and M70, respectively (See also 4^ in Table [TJ. The 
spatial size of each box is —5 X 10® cm < x' < 5 X lO^cm and 
< 2 < 1 X lOl'^cm. The upper and lower panels correspond to 
M150, Mref, M70 and M50 models. Note that some panels in M150 
and Mref are lacking since we stop the simulation at the time of 
breakout. 

identical among all models and follows a rather spheri- 
cal contraction, since the centrifugal force plays a minor 
role. Note that we find that the density distribution at 
the inner boundary is slightly oblate but this does not 
affect the subsequent jet propagation although it might 
have consequences for the jet production (which is not 
properly simulated here). 

During the jet propagation phase, on the other hand, 
the results of jet evolution are very different among each 
model (see Table [T] and Figure [2]). We also show that 
the summary of our results in Tabled] For model M50, 
the forward shock wave does not move out and almost 
stagnates around the inner boundary despite the success- 
ful operation of the central engine {ti in Table [1] denotes 
the time of initiation of central engine). In fact, no colli- 
mated feature can be seen for M50 in the lowest panels in 
FigurcO This is attributed to the fact that the jet power 
does not exceed the ram pressure of the inflowing mate- 
rial, and the forward shock wave stagnates or is advected 
inwards. For models with successful jet breakout, on the 
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Fig. 3. — From the top to the bottom, hemispherical neutrino lu- 
minosity, mass accretion rate, black hole mass, spin parameter and 
conversion efficiency from accretion energy to neutrino luminosity 
(77 = Lj/Mc^) as a function of time from the onset of the col- 
lapse. Each color represents different models. Note that the model 
name with "extend" indicates the extrapolate results by analytic 
treatment (see text for more details.). 



increases with time (see the 4th panel in Figure [3]), which 
produces a powerful jet as a result of the large neutrino 
deposition energy rates. It is also important to note that 
the onset timing of central engine (ti) also greatly affects 
the outcome of explosion. As shown in Table [1] and il- 
lustrated in Figure [2J the jet production is significantly 
delayed for a slower rotational model. This is due to 
the neutrino luminosity being weaker for both smaller 
accretion rates and larger black hole masses (see Eq. ^ 
for the dependence of m and M^h)- In fact, for M50, 
although the Kerr parameter reaches > 0.9 at the end 
of our simulation, the neutrino luminosity is not large 
enough to move out the forward shock wave. 

According to these results, we infer that neutrino 
driven jets may not penetrate progenitors with extended 
envelopes, since significant large mass might be able to 
accrete to the black hole before jet breakout. The weak 
neutrino luminosity resulting for a sizable increase in 
black hole mass could result in the jet becoming non- 
relativistic ejecta. Therefore a compact progenitor is an 
inevitable requirement for a successful of neutrino-driven 
jet breakout. For the massive envelope progenitors such 
as PoP III or Red (Blue) supergiants, a different process 
m ight be required to p owered GRBs (see also discussions 
in iSuwa fc lokal (|20Tlh '). 

The foremost important result in this study is that the 
jet succeeds to break out from the star except for M50, 
which corresponds to the model with the slowest rotation 
rate. The neutrino-driven jet with 6op ~ 9° produced in a 
rapidly rotating compact Wolf-Rayct star can potentially 
give rise to a GRBs. We also find that the time-averaged 
accretion-to-jet conversion efficiency among successful jet 
breakout models is roughly 77 ~ 10"'^, while ry for M50 
can not reach 10~^ and never accomplish the jet break- 
out. This result is roughly consistent with our previous 
work jNagakura et al. 2012). The analytical criteria in 
(jNagakura et al.il2012D also shows the opening angle de- 
pendence for the successful jet breakout, which is the 
threshold 77 increases with (^op)^- Therefore, according 
to this result, we would like to give an important caution 
that jets with wider opening angle are more difficult to 
penetrate the star than the present results, which indi- 
cates that the threshold progenitor rotation rate differs 
in accordance with the jet opening angle. 



other hand, the jet also can not move forward quickly af- 
ter the initiation of the central engine. However, due to 
an increase in the Kerr parameter of the black hole over 
time, the jet power eventually exceeds the ram pressure 
of the inflowing material (Figure Once the forward 
shock wave is able to move out, the jet interacts with 
the stellar mantle and gives rise to a cocoon. The hot 
cocoon helps jet confinement and helps to preserve the 
jet's strong outgoing momentum and energy flux, even- 
tually leading to a successful breakout. 

Figure [3] shows the evolution of hemispherical neutrino 
luminosity, mass accretion rate, black hole mass, Kerr 
parameter and conversion efficiency from accretion en- 
ergy to neutrino luminosity {t] = Lj/Mc^) as a function 
of time from the onset of the collapse. The chief cause 
for the different jet propagation behavior is the sensitive 
dependence of the neutrino luminosity on the Kerr pa- 
rameter (see Eq. For the fast rotational model, the 
angular momentum of the black hole is very large and it 



3.2. Post breakout phase: Analytical Formula 

As we have already mentioned, our numerical simula- 
tions are terminated at the time of jet breakout. How- 
ever, it is interesting to extend the result of our numerical 
simulations to post break out stage, which allows us to 
estimate the expected observational differences among 
the computed models (see Section 13. 3p . We employee 
the following analytic approximations in order to see the 
subsequent evolution after the jet break out; 



"''■(0) 

r dM 

Mbh^Mhhib)+ —^{r)dr (4) 

"''■(6) 
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where 



Mffir) 



GMr 



dMr/dr 
dtff/dr 
1 SirGMrtffp 
'iP 3Mr - inr^ p' 



(5) 



(6) 



The above analytical estimati on is essentially simila r 
to the approach p resented in iNagakura et al.l ()2012D : 
iSuwa fc lokal (|2011[ ). The free fall time tff{r) can be 
determined by based on the assumption of spherical sym- 
metric envelope contraction. Here, th denotes the time 
at jet break out. Note also that the functions of Mr{r) 
and p{r) are extracted from the table of 16TI model, r?, 
is determined by the assumption of Mbh{tb) = Ajfr(?'(b))- 
The non-dimensional parameter /3 is determined to en- 
sure that M{tb) is equal to Mff{r{p)). According to this 
procedure, the neutrino luminosity and mass accretion 
rate can be smoothly connected from the results of nu- 
merical simulations. The evolution of the Kerr param- 
eter is determined by Eq. ([1]). Note that, since model 
M50 does not succeed the jet breakout, there is no post 
breakout phase for this model. Results of these analytical 
extension are also described in Figure [3] Figure S] shows 
the comparison between the results of extended numer- 
ical simulations and analytical estimation for the refer- 
ence model. The extended numerical simulations are per- 
formed until t = 40s, which corresponds to about ten sec- 
onds after the jet breakout. Note that we do not broaden 
the computational region, since the stellar contraction is 
not affected by the jet dynamics in the outer parts of 
the star. As shown in this figure, the time evolution of 
black hole mass {Mbh), Kerr parameter (a) and conver- 
sion efficiency (ry) are almost identical between results 
of the extended numerical simulation and the analytical 
calculation. For luminosity (Lj ) and mass accretion rate 

(A/), on the other hand, the analytical calculations are 
slightly larger than the results of our numerical simula- 
tions. This may be attributed to the fact that the ana- 
lytical approach neglects stellar rotation, which increases 
the mass accretion rate, and consequently overestimates 
the neutrino luminosity. However, these differences are 
within ten percent. Therefore, we confirm that the above 
analytical approach qualitatively well describes the time 
evolution of the dynamics of the jet. In the following sub- 
section, we discuss the observational consequence with 
the aid of the analytical approach in the post breakout 
phase. 

3.3. Observational consequences 

We first divide the energetics of the jet into two parts, 
which are relativistic jet component (Ej) and the co- 
coon c omponent {Edg) (see lSuwa fc lokal (|2nilh : lMatzned 
(|2003l )). The Ej is calculated based on the assump- 
tion that all the injected energy after the jet break out 
goes into the relativistic component, i.e, Ej is given by 
j^{Lj/2)dt. Note that 1/2 factor comes from the as- 
sumption of equatorial symmetry. 

As shown in Table [1] Ej increases with an increase of 
stellar rotation. For the purpose of studying the outcome 
of explosion in more detail, we further divide the energy 




0'-0001 



Fig. 4. — Same as FigureO but comparison between results from 
extended numerical simulations and analytical approaches for the 
reference model. 

of relativistic jet into Ej^^^g and Ej^^^g ^. Ejy^^^ is cal- 
culated by the same manner as Ej except for the integra- 
tion is carried out when Lj/2 > 10^°erg/s, while Ej^L^^ ^ 
is calculated with a condition Lj/2 > 5 x 10'*^erg/s. In 
addition, we also list the corresponding time duration 
of each component as Tjy^^g and TjyL^^^ in Table [TJ 
In all cases, these timescales are several tens of seconds, 
which are comparable with typical time scale of prompt 
phase of GRBs. It should be noted, however, according 
to Br omberg et al.. ^OlJ), the duration of prompt phase 
of GRBs may be modified by the durati on of jet penetra- 
tion, i .e, tj = te—tb (sec also Eq. (2) in lBromberg et aTl 
()2012|) ). where t^, te, h denote observed duration of 
the prompt phase, the central engine working time and 
the duration of the jet penetration phase, respectively. 
Therefore, the actual observed duration will be smaller 
than the Tj, and it is substantially modified especially 
for slower rotation models (e.g. M70) (since % is larger 
with slower rotation). It is also interesting to note that 
model M70, which is the slowest model among successful 
jet breakout models, produces the weak explosion and 
does not have i?j>L5o due to the low luminosity of the 
jet. Note also that this luminosity is upper limit for ob- 
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Hypernovae. 

We further calculate the isotropic energy {Ej(^iso)) for 
Ej, EjyL^g^ and Ej^L^g, and also isotropic peak lumi- 
nosity . which are shown in Figure [5] In this cal- 
culation, the jet opening angle is assumed to be f 
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which is the same as the root of injected jet in our sim- 
ulations. Note again that we neglect the conversion ef- 
ficiency from hydrodynamic energy to radiation, so our 
results are still at the qualitative level and give the only 
upper limit of GRB radiation. In addition, the time 
evolution of neutrino luminosity does not capture the 
rapid time variability of central engine since our numer- 
ical simulations do not incorporate the black hole ac- 
cretion disk system. According to these ambiguities, 
the time evolution of neutrino luminosity as described 
in the upper panel of Figure [3] are different from the ob- 
served light curve in reality. It should be noted, however, 
that our analysis is meaningful to give the constraint 
the neutrino luminosity and energy as the upper limit. 
Based on the above assumption, we find that -E'j(iso) is 
~ lO^^erg, while Lp(^iso) is ~ lO^^erg/s, which are suffi- 
ciently large to explain GRBs. We would like to point 
out that, for the jet with for rapidly rotating progenitor 
(M150), large fraction of energy are radiated with high 
luminosity jet {Lj > 10^°erg/s), while more than half of 
jet energy for M70 would be radiated with low luminous 
jet (Lj < 5 X 10*^erg/s). According to these results, we 
suggest that neutrino-driven jet is capable of producing 
several types of bursts by the different progenitor rota- 
tion, which may be the origin of observational different 
bursts such as GRBs, XRFs. The failed GRBs would be 
also explained by neutrino-driven central engine when 
the progenitor is slowly rotating. 

4. SUMMARY 

We present the numerical results of neutrino-driven jet 
propagation in a rotating Wolf-Rayet star. By changing 
the rate of progenitor rotation, we discuss jet penetrabil- 
ity and their observational consequences with the aide of 
analytic extrapolation in the post breakout phase. We 
show that every model except for M50 succeeds to break 
out the star. Especially, Mref and M150, which corre- 
spond models with sufHcient rapidly rotation, have the 
relativistic outflow component as ipj-j^o) ^ lO^^erg/s and 
^j{iso) ^ lO^^crg, which are sufhcicntly large to explain 
GRBs. On the other hand, the energy in the cocoon com- 
ponent Edg is ~ lO^^erg for models with successful jet 
break outs, although it remains an open question whether 
the jet or the cocoon expansion could give rise to enough 
nickel production to explain the GRBs-Hypernovae con- 
nection. One of the other important results in this study 
is that model M50, which corresponds the model with 
slowly rotational model, can not succeed the jet break- 
out (failed GRB). Therefore, there is the threshold SAM 
distribution between model M70 and M50 for the success 
of jet penetration. It should be noted, however, that the 
threshold rotation, no doubt, strongly depend on the jet 

opening angle, and our results are only adequate for the 

sions of cocoon propagatio n in iMeszaros fc ReesI (|2001h : canonical jet opening angle, 6op = 9°. Although there are 



Fig. 5 . Ej^2so) tind Lp(^iso) relation among models with the suc- 
cessful jot breakout. -Ejjiao) is calculated from Ej (Red), Ej^^^^ 
(Green) and Ejyi^^^ (Blue), respectively. 

served luminosity since we neglect the conversion process 
from hydrodynamical energy to gamma-rays. According 
to these facts, we infer that the neutrino-driven jet from 
the compact Wolf-Rayet star, whose rotation rate is be- 
tween model M70 and M50, may produce very low lumi- 
nous type burst, possibly LLGRBs. These results are 
qual itatively consistent w i th th e previous stud i es (se e 
e.g. iLopez-Camara et"al] ()2009[ ): iLindner et all (|2010( ): 
iMilosavlievic et al.l ( 2012[ )). It should be noted, however, 
that some LLGRBs are observed with the extreme long 
duration, which populations can not be explained by the 
results of current studies. Therefore, the neutrino-driven 
jet may contribute to only LLGRBs with typical dura- 
tion of prompt burst ('^ 10s). 

On the other hand, the energy of cocoon component 
can be estimated as th e diagnostic energy at t = tt (see 
iNagakura et al.l ( 2012|) ). It is important to note that, as 
shown in Table [11 Edg is typically ~ lO^^erg for mod- 
els with successful jet breakout. This may be attributed 
the fact that the jet with a slowly rotating core tends 
to be weaker power and it spends longer time to pen- 
etrate the star. Therefore, despite its low jet luminos- 
ity, a large fraction of jet energy has been consumed for 
sweeping aside the stellar mantle and accumulated as the 
cocoon energy, which eventually reaches ~ lO^^erg. The 
cocoon material is expected to contribut e the subsequent 
explo s ive event after the p rompt phase (jTominaga et al] 
120071: iLazzati et al] I2012D and at the afterglow phase 
( Ramirez-Ruiz et al The results of neutrino- 

deposition presented in this paper are not able to dis- 
cern whether sufficient large nickel production might take 



place to explain hypernova ex plosions (see also discus 

cocoon propagatio n in IMeszaros fc ReesI (|2001h : 

iRamirez-Ruiz et al.l (|2002D : iMatzneij (|2003[) ). If nickel 



is not effectively produce at the jet interaction region 
alternative pathways such as the disk wind by viscous- 
heating or magnetic-driven wind from central engine 
would be required to explain the link between the GRBs- 



some limitations in this study, we suggest that neutrino- 
driven jet is capable of producing several types of bursts 
(and also include failure branch of burst i.e, failed GRBs) 
by the different progenitor rotation. 
Finally, we would like to note that the results presented 
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in this paper are optimistic. As one of the large uncer- 
tainties, the actual mass accretion rate would be smaller 
than the results obtained in this paper, since some frac- 
tion of the mass is expected to escape from the disk 
rather than being accreted to th e black hole due to neu- 
trino winds or viscous heating (jMacFadven fc WooslevI 
I1999D . In addition, the increase rate of Kerr parameter 
would be slower than the current result, since the SAM 
of the inflow matter is smaller than the ISCO due to the 
effect of pressure gradient. Note also that the disk wind 
also extract the angular momentum of accretion matter. 
Since the neutrino deposition rate depends sensitively on 
the mass accretion rate and Kerr parameter, the jet dy- 
namics would be affected by these effects. Note also that 
the neutrino-driven jet can not explain the extreme long 
duration of bursts and other populations are necessary 



to explain these peculiar events. The other factors such 
as viewing angle may also cause the obs ervational differ- 
ence among GRB population (See e.g. lYamazaki et al.l 
(|2002l 120031) : iGranot et all (|2005D ). More quantitative 
discussions will be conducted in our forthcoming paper. 
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